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Abstract. - Rotational states for trapped bosons in an optical lattice are studied in the framework 
of the Hubbard model. Critical frequencies are calculated and the main parameter regimes are 
identified. Transitions are observed from edge superfluids to vortex lattices with Mott insulating 
cores, and subsequently to lattices of interstitial vortices. The former transition coincides with 
the Mott transition. Changes in symmetry of the vortex lattices are observed as a function of 
lattice depth. Predictions for experimental signatures are presented. 
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Ultracold atoms in optical lattices combined with trap- 
^ping potentials offer rich possibilities for exploring new 
i physics. In particular, the presence of a trap allows for 

the coexistence of phases such as the Mott insulator and 
' sjuperfluid, as was predicted by Jaksch et al. [1], realized 
f—^by Greiner et al. [2], and directly verified by two inde- 
0\ pendent groups in 2006 [3,4]. The harmonic confinement 

results in a shell-like structure, where superfluid and Mott 
C^l insulating regions occupy alternating intervals (in one di- 
i^^mension), annuli (in two dimensions), or spherical shells 

(in three dimensions). 
OO , 

, Motivated by the progress in the study of vortex physics 
L^in trapped condensates, rotational states in optical lattices 
• — Ijave attracted some recent attention. At one end of the 
rich spectrum of possibilities is the case of a weak pcri- 
?H qdic potential that acts as a lattice of pinning centers for 



vortices; this has been realized in experiment using a coro- 
tating periodic potential [5,6]. For deeper optical lattices, 
one may consider a discretized description of the gas, the 
bosonic Hubbard model [1,7], and study how rotational 
states interact with the quantum phases of the model. It 
has been predicted that quantized vortices in the Hubbard 
model possess a Mott insulating core [8,9], and also that 
in the limit of a deep lattice, rotational states appear as 
an edge superfluid surrounding a Mott insulator [10]. In 
the limit of rapid rotation, when the number of vortices is 
comparable to the number of lattice sites, intriguing phe- 
nomena are predicted such as strongly correlated vortex 
liquids [11] and fractal densities of states, the so-called 
Hofstadter butterfly [9,12-14]. 

This study aims to investigate the interplay between 
the quantum phase transition, rotation, and trapping for 



atoms in a optical lattice described by the bosonic Hub- 
bard model, and to describe the different types of rota- 
tional state that can be found in such a system. It will 
be seen that at moderate rotational velocities, the rota- 
tional states of the system are closely connected to the 
Mott transition. Since much of the qualitative physics is 
expected to be brought out already in two dimensions, the 
study is carried out in 2D. 

A gas of spinlcss bosonic atoms in an optical lattice is 
known to be well described by the Hubbard model [1] . In 
the presence of rotation, the Hamiltonian is [8,10] 
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where i,j label the lattice sites, < ij > denotes near- 
est neighbors, and rj is the position vector of site i 
with length r^. U is the on-site interaction strength, 
t is the tunneling matrix clement, and /i is the chemi- 
cal potential. The rotation introduces a phase <f>ij. An 
alternative Hamiltonian was proposed in Ref. [7] using 
a corotating ansatz for the single-particle orbitals; how- 
ever, the quantitative differences would be exceedingly 
small in the present case. We state energies in units 
of the recoil energy En and distances in units of the 
lattice constant d = A/ 2 where A is the wavelength of 
the lattice beams, and for definitcness the Wannicr func- 
tions are approximated by Gaussians. In this approxi- 
mation t = Vo 3/4 7r 5 / 4 [l - (2/tt) 2 ] cxp(-7r 3 / 2 VVb/4)/4, and 



U = Vr 



3/4, 



./\/2, where Vq is the depth of the lattice po- 
tential and a is the s-wave scattering length of the atoms. 
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For the calculations, we have assumed Rb atoms with a 
scattering length a = 5.77nm and a square lattice made 
out of two pairs of counterpropagating beams with wave- 
length A = 780nm. We define ft = TtKl/^/2, where f2 is 
the rotational frequency, and u> — Trhuj/^/2, where ui is 
the frequency of the harmonic magnetic trapping poten- 
tial. Finally, the phase (f)^ = (m/h) f* J dr-(ilxr), where 
in this 2D study the vector f2 points out of the plane. 

The Hamiltonian in Eq. (1) describes an optical lattice 
that rotates with the frequency f2. Such a potential was 
realized in Ref. [6]. Other ideas for creating a rotating 
state in an optical lattice include exciting the desired vor- 
ticity in a trapped condensate, and applying the optical 
lattice afterwards. Alternatively, one may stir the gas by 
using a timc-dependent optical or magnetic potential [15], 
or by phase engineering [16-18], in the presence of a static 
optical lattice potential. However, the latter options are, 
rigorously speaking, not described by the Hamiltonian in 
Eq. (1). 

We now map out the rotational phases of the trapped 
system. The large number of parameters makes for a rich 
phase diagram. The rotation with frequency f2 gives rise 
to a centrifugal potential that competes with the confining 
potential. At the critical value ft = lj, the centrifugal po- 
tential balances the trap potential, and thus this criterion 
gives the upper bound for the angular velocity. However, 
the present method of calculation is not reliable in the 
fast rotating limit, so we will only consider the regime 
well away from this upper bound. In addition, we work 
with a rather weak trapping potential, Co = 0.025, and 
hence the scaled rotational frequency O must be equally 
weak. Invoking the analogy with a system in a magnetic 
field, (l is the flux through each plaquette [19]. Thus, the 
present study takes place in the limit of weak flux, preclud- 
ing effects such as Hofstadter-butterfly physics [9,12-14]. 
However, this limit exhibits a number of interesting phys- 
ical effects connected with the Mott transition, which will 
be the subject of the present paper. Furthermore, we deal 
only with moderate densities, keeping the chemical poten- 
tial fi = 0.5U for the most part and increasing it to maxi- 
mally fi = 2.0U. In this way we can focus on vortex-lattice 
physics without invoking unnecessary complications such 
as alternating Mott and superfiuid shells. Thus, at fixed Co, 
we now explore the parameter space spanned by moderate 
ranges of Vq, fi, and fi. 

In the absence of rotation and external potential, the 
Hamiltonian (1) undergoes a quantum phase transition be- 
tween two types of ground state at a critical value of the 
ratio of the tunneling and interaction parameters, t/U, de- 
pending on the chemical potential fi [20]. When the ratio 
t/U is large enough, there is phase coherence over the en- 
tire sample which puts the system in its superfiuid state, 
and a Bose-Einstein condensate is formed. For weaker 
tunneling, phase coherence is lost, number fluctuations are 
suppressed, and the number of atoms per site is locked to 
an integer. This is the Mott insulating state. The critical 
tunneling in two dimensions at unit filling was in Ref. [21] 



calculated to be (t/U) c = 0.060, while in the mean field 
approximation, which we shall use here, one finds the lower 
value (t/U) c = 0.042. In a trapping potential V(r), su- 
perfluids and Mott insulators coexist in spatially separated 
regions, where the local phase is determined by the local 
chemical potential /i; oc = fx — V(r) [1,3,22]. The size of 
the regions depends on the ratio t/U, such that for large 
enough t/U, the whole sample is superfiuid, and as t/U 
decreases, the Mott insulating regions grow. 

If the atoms are set into rotation, it is the superfiuid 
fraction that rotates since the Mott fraction is insulating. 
A continuous superfiuid, or a trapped Bose-Einstein con- 
densate in the absence of an optical lattice, responds to 
rotation by forming a lattice of quantized vortices with 
hollow cores [23]. There exist excited states confined to 
the cores, but in a weakly interacting sample their pop- 
ulation is modest [24]. In a numerical study of an ex- 
tended Hubbard model, with interactions between nearest 
neighbors, it was found that the superfiuid forms vortices, 
and the cores are Mott insulating [8]. In other parameter 
regimes, interstitial vortices may form, with only modest 
suppression of the condensate density [9,11,13]. Yet an- 
other parameter regime was investigated in Ref. [10]. If 
almost all of the sample is Mott insulating, rotation will 
result in an edge superfiuid containing only a few rotating 
atoms, and the bulk is unaffected. 

Numerical calculations have been carried out on a 
100x100 lattice in the mean-field Gutzwiller approxima- 
tion, which has proven to be reliable in a range of situa- 
tions [1,25,26]. To solve for the rotating states, a phase 
winding has been imposed onto the nonrotating solution 
whereafter the energy has been minimized using a con- 
jugate gradient method [10]. 5 single-particle states per 
lattice site were allowed for in the calculations. Solutions 
for the case of two phase singularities, i.e., a total phase 
winding of Air around the circumference of the sample, are 
shown in Fig. 1. The Gutzwiller approximation allows for 
direct calculation of the condensate wavefunction Zi — (a/) 
at each point; its square na = \zi\ 2 is denoted the di- 
mcnsionlcss condensate density and is plotted alongside 
the dimcnsionlcss density, which is defined as the mean 
number of atoms on the site, 7ij = (ajaj). The chemical 
potential is chosen as fi = 0.5U. The mean-field critical 
point for the Mott transition at t/U = 0.042 corresponds 
to Vo = 27.2. At Vq = 26.2, the entire system is in the su- 
perfiuid state, and the condensate density is nonzero but 
contains two depressions which coincide with phase sin- 
gularities; these are vortices. At these points, the total 
density is unity, supporting the conclusion of Ref. [8] that 
the cores are filled with Mott insulating atoms. As Vo is 
increased, the condensate density in the center decreases 
and the edge-superfluid state is entered. At the Mott tran- 
sition point Vq = 27.2, the condensate density in the center 
is nco = 3.2- 10 -6 , and the two phase singularities are sep- 
arated by a single lattice site. For larger Vq , «co vanishes 
to within numerical precision close to the center, so that 
the phase varies erratically. Clearly, the Mott transition 
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Fig. 1: Density and phase profiles for a system containing two 
phase singularities, as the Mott transition is crossed, (a), (d), 
(g), (j): Cross-section at y = of the dimensionless density 
rii (upper curves) and condensate density nci (lower curves), 
(b), (e), (h), (k): Condensate density. Bright shades signify 
high density and vice versa, (c), (f), (i), (1): Phase of the 
condensate wavefunction Zi, increasing with brightness from 
(black) to 2-7T (white). The lattice potential height is chosen as 
V = 27.2 (a-c); V = 26.7 (d-f); Vo = 26.2 (g-i); and V = 22 
(j-1). This corresponds to t/U = 0.0418, 0.0447, 0.0478, and 
0.0868, respectively. The chemical potential is fi — 0.5C7 and 
the trapping frequency is u> = 0.025. The rotating frequency is 
n = 0A7Q in (a-i) and fi = 0.63d; in (j-1). 



is accompanied by a transition between the two types of 
rotational state. 

Since the Mott transition is second order, there is no 
surface tension associated with the phase boundary and 
hence no tendency for the vortices to stick together and 
form giant vortices. In the absence of a lattice poten- 
tial, a giant vortex is never stable in a harmonic trap [27], 
and apparently, the same holds for the trapped Hubbard 
model, as long as the whole system is in the supcrfluid 
phase. 

Moving further from the Mott transition, i.e., decreasing 
the lattice depth, the vortex cores are observed to decrease 
in size. When a vortex core is smaller than the lattice con- 
stant, one enters the regime of interstitial vortices, where 
the depletion of the condensate is modest. An example 
is shown in Fig. 1 - The reason is that this system is 
discrete, and the vortices naturally tend to form in the in- 
terstices, just as observed in most lattice models. In other 
words, on a plaquette of 2 x 2 lattice sites the phase may 
wind by 2ir, and the condensate density does not have to 
vanish on the sites because the phase singularity is in the 
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Fig. 2: (a) Size f e fj of a vortex core as the lattice depth Vo, 
and hence the tunneling over on-site energy, t/U, is varied, 
(b) Ratio of core to bulk condensate density, nco/ncb, as a 
function of Vo and t/U. The chemical potential is fi = 0.5(7, 
the trapping frequency u> — 0.025, and the rotational frequency 
is fl = 0.5(D. 



void between them. A discernible vortex core forms when 
the balance between the interaction energy and the cen- 
trifugal energy is such that a large vortex core is favored. 
In a trapped Bosc-Einstcin condensate, the size of a vor- 
tex core is given by the healing length, £ = y/t/(Uncb), 
where neb is the mean condensate density away from the 
vortex [23]. In the vicinity of the Mott transition, the 
formula cannot be trusted quantitatively, but the depen- 
dence on neb provides a qualitative explanation of why 
the vortex core size decreases when one moves away from 
the Mott transition. 

The size of the core of a single vortex has been calculated 
by fitting the central part of the condensate density profile 
nci to the the interpolation formula [23, 28] 
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where i = (i x ,i y ) is the Cartesian coordinate of a lattice 
site and i x o, i y o is the location of the center of the vortex. 
Two series of calculations are made, one where the chem- 
ical potential is fixed at [i = 0.5/7 and the lattice depth is 
varied; and one where fi is varied with the lattice depth 
fixed at Vq = 25. The results for nco and £ e // as functions 
of t/ U is shown in Fig. 2, and as functions of fx in Fig. 3. 

For simplicity, the angular frequency has been fixed at 
Q = 0.5a), although the central singly quantized vortex is 
not the global energy minimum over the whole range. 

The crossover from Mott-core vortices to interstitial vor- 
tices is most clearly seen in Fig. 2(a), where the conden- 
sate density in the vortex core nco is recorded. In the 
Mott phase i.e., for t/U < 0.042, the condensate den- 
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Fig. 3: (a) Size £ c ff of a vortex core as a function of chem- 
ical potential /j,. (b) Ratio of core to bulk condensate den- 
sity, nco/ncb, as a function of p. The lattice depth is kept 
at Vo — 25Er, which implies t/U = 0.0564, the trapping fre- 
quency Q = 0.025, and the rotational frequency is Q = 0.5a). 



sity in the center is zero to within numerical accuracy. 
As soon as the superfluid phase is entered, nco increases 
to a small but finite value. However, there is a second, 
reasonably sharp transition, where nco abruptly rises to 
a value of about half the condensate density outside the 
core. This transition is seen to occur when the fit pa- 
rameter £ e ff has dropped to a value well below one lattice 
site. The abrupt transition suggests that we can regard 
the interstitial vortices and Mott-core vortices as distinct 
rotational phases. In Ref. [9], a homogeneous system with 
periodic boundary conditions and a large vortex density 
was studied. A similar abrupt transition was found there 
between vortices centered on sites and vortices centered 
on plaquettes. However, as was noted in Ref. [9], since 
energy differences between different types of ground state 
may be small close to the Mott transition, it may very well 
be that the exact phase diagram is substantially different 
from that presented here. 

In Fig. 3, the slightly oscillatory dependency of nco as 
a function of /j, reflects the proximity in phase space to an 
alternating sequence of Mott insulating regions centered 
around /j, = 0.5,1.5, etc. [20,22]. It is clearly seen that 
away from the Mott transition, the core size £ e ff decreases 
and the condensate density in the core nco jumps to a 
value of the same order as the bulk condensate density, 
but close to the Mott insulating regions it is significantly 
suppressed. 

As the angular velocity is increased, the minimum- 
energy state of the system contains an increasing number 
of vortices. In Fig. 4 minimum-energy configurations with 
six phase singularities are plotted far from and close to 
the Mott transition. In the superfluid phase, a rudimen- 
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Fig. 4: Density and phase profiles for a system with six phase 
singularities. The quantities displayed are as in Fig. 1. (a- 
c): V = 27.4, Q,/u = 0.68, t/U = 0.0407. (d-f): Vo = 27.2, 
fi/w = 0.68, t/U = 0.0418. (g-i): Vo = 26.7, fi/w = 0.70, 
t/U = 0.0447. (j-1): Vo = 20, Sl/u) = 0.77, t/U = 0.1176. The 
chemical potential is /j, — 0.5(7. Note that the parameters are 
slightly different from those in Fig. 1. 

tary vortex lattice is formed, resembling those found in 
trapped condensates (as well as superconductors and liq- 
uid helium, of course). The vortices are filled with Mott 
insulating atoms so that the total density profile is not 
strongly affected. It is seen that the lattice changes from 
fivefold to sixfold symmetry for Vo > 22; this is further 
discussed below. At the Mott transition, Vq = 27.2, the 
phase singularities do not seem to have molten together 
- the shortest distance is about 5 lattice sites - but the 
central density is nco = 8 ■ 10 -6 , and as Vo is further in- 
creased, the phase begins to vary erratically in the center, 
as seen in Fig. 4(c). The figure supports the conclusion 
that the Mott transition marks the transition between a 
vortex-lattice state and an edge superfluid. The rotating 
force does not seem to shift the critical point at least for 
the rather modest values of O considered here. 

The critical frequencies Qc,q for the thermodynamic sta- 
bility of a state containing q phase singularities is cal- 
culated by comparing the energies of different stationary 
solutions of the Hamiltonian (1), computed by choosing 
different phase windings as the initial condition for the 
energy minimization. The result for the lowest values of q 
is shown in Fig. 5. It is known from trapped Bose-Einstcin 
condensates that the critical frequency does not directly 
decide whether a rotating state can be created, since that 
creation depends on stability properties of surface exci- 
tations [29]. For a complete picture, the phase diagram 
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Fig. 5: Critical frequencies 0.c, q for rotational states in the 
Hubbard model, as a function of the optical lattice potential 
depth Vo. The chemical potential and external trapping po- 
tential are fixed at fi = 0.5(7 and ui = 0.025. The numbers 
denote the total circulation of the state that is thermodynami- 
cally favorable in the indicated interval. The dotted line is the 
prediction of Eq. (3). The horizontal dashed lines mark the 
transitions between the rotational phases (see text). The dots 
mark the transitions between vortex lattices of different sym- 
metry, for a total of 6, 7, and 8 quanta, respectively, in order 
of increasing Vo- 



shown in Fig. 5 should be supplemented with an analysis 
of the dynamical process for exciting vortices, but such an 
analysis is not attempted here. 

Some features of the phase diagram can be understood 
quantitatively. Deep in the Mott phase, rotation results 
in an edge current. Since the supcrfluid forms an annulus 
around the central Mott insulating core, the energy cost 
for creating a circulating current - a phase winding - is 
small and the critical angular velocity for exciting such a 
current will be accordingly small. In order to estimate the 
critical angular velocity, assume that the rotation does 
not appreciably change the number of superfluid atoms, 
but only induces a phase winding. The superfluid atoms, 
whose number we denote Nc, reside in a thin shell of some 
thickness d whose radius is R, the size of the system. This 
radius is obtained from the condition that the external 
potential matches the chemical potential: in dimensionful 
units, R 2 = 2fx/[m(tu 2 — fi 2 )]. The velocity at radius R 
is v — (h/m)q/R, where q is the phase winding. The 
energy associated with the rotation is just equal to half 
the squared velocity times the number of superfluid atoms, 
E TO t = Ncmv 2 /2, while the total angular momentum is 
L = NcmvR. Hence the angular velocity for stabilization 
of q units of angular momentum is given by the solution 



to the equation 
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where we have returned to dimcnsionless units. It is seen 
in Fig. 5 that the estimate holds very well for the first 
critical frequency. However, it should be noted that the 
very close quantitative agreement is probably fortuitous, 
since the higher critical frequencies do not agree to the 
same level of precision. 

For a Bose-Einstein condensate, it is known that all the 
critical frequencies should approach the trap frequency 
as the ratio of kinetic energy to interaction energy is 
increased [23]. The reason why the critical frequency 
docs not increase as Vq is decreased below approximately 
Vo = 20 is that the system enters the interstitial vortex 
phase, where estimates for Bose-Einstein condensates can- 
not be trusted. Besides, this study is performed at a fixed 
ratio fi/U, and the total number of bosons increases as 
a function of Vo at fixed f2. The transition between the 
Mott-core and interstitial vortex phases at Vo = 22.5 is 
indicated in Fig. 5. 

The solid circles in the phase diagram mark the tran- 
sitions from ring configurations of vortices, for large Vo, 
to lattice-like vortex arrays for smaller Vo, as exemplified 
in Fig. 4. This transition is most likely due to competi- 
tion between two energies. The tendency for the vortices 
to arrange themselves on a ring arises from the "acciden- 
tal" fact that at the chosen chemical potential fj, = 0.5, 
the condensate density nc develops a local minimum at 
a finite radius, which attracts the vortices. Further from 
the Mott transition the minimum weakens, and in combi- 
nation with the repulsion between vortices it leads to the 
more familiar vortex lattices. By adding rotation to the 
already complicated phase structure of trapped bosons in 
optical lattices, a wealth of such different vortex configu- 
rations can be expected. 

The detection of these intriguing rotating states was 
discussed in Rcf. [10], where it was observed that time- 
of-flight imaging can yield quantitative information about 
the state. The image of the atom cloud following a re- 
lease from the trap will be distinctive: the Mott insulat- 
ing atoms expand incoherently, but the Bose-Einstein con- 
densed fraction assumes a shape close to the Fourier trans- 
form of the condensate wave function Zi\ it is seen that 
the Fourier transform contains plenty of information on 
the original spatial distribution. Because of the periodic 
potential, copies of this image appear translated by recip- 
rocal lattice vectors [2] . Six examples of expected time-of- 
flight images are shown in Fig. 6. In order to enhance con- 
trast, the time-of- flight technique may be combined with 
selective removal of Mott insulating atoms [3]. Time-of- 
flight imaging can thus easily discriminate between vortex 
lattices and edge superfluids. It will also directly reveal the 
symmetry of the vortex lattice. The transition between in- 
terstitial and Mott-filled vortices cannot, however, be eas- 
ily seen in this way. in order to spot this transition, the 
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Fig. 6: Examples of density profiles after expansion of a rotat- 
ing cloud, (a) is the expected time-of-night image of the state 
corresponding to Fig. l(a-c); (b) corresponds to Fig. l(g-i); (c) 
corresponds to Fig. l(j-l); (d) corresponds to Fig. 4(d-f); (e) 
corresponds to Fig. 4(g-i); and (f) corresponds to Fig. 4(j-l). 



method of direct detection of Mott insulating sites [3,4] 
should be used instead. 

In conclusion, the rotating states for trapped bosons 
in an optical lattice have been mapped out. Three basic 
types of state have been identified, namely edge current, 
Mott- filled vortices and interstitial vortices. The transi- 
tion between the two former phases coincides with the 
Mott transition. The transitions between the three types 
of state are found to be quite abrupt. Critical frequencies 
are calculated for low- lying rotating states. It is likely that 
both time-of-flight experiments and occupation-number 
sensitive detection are needed in order to detect the states 
discussed here; examples of predicted experimental sig- 
natures are calculated. Seeing that symmetries of vortex 
arrays are sensitive to the radial dependence of the con- 
densate density, it is conceivable to think of vortex pat- 
terns as an experimental probe for the phase diagram of 
the trapped bosonic Hubbard model. 

This research was conducted using the resources of High 
Performance Computing Center North (HPC2N), and fi- 
nancially supported by the Swedish Research Council, 
Vetenskapsradet. 
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